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Abstract. The metallicity evolution and ionization history of the universe must leave its imprint on the Cosmic Microwave 
Background through resonant scattering of CMB photons by atoms, ions and molecules. These transitions partially erase origi- 
^\ ' nal temperature anisotropies of the CMB, and also generate new fluctuations. In this paper we propose a method to determine 

the abundance of these heavy species in low density (over-densities less than lO"* - 10^^) optically thin regions of the universe 
CN ■ by using the unprecedented sensitivity of current and future CMB experiments. In particular, we focus our analysis on the 

^ I sensitivity of the PLANCK HFI detectors in four spectral bands. We also present results for 1=220 and 810 which are of interest 

for balloon and ground-based instruments, like ACT, APEX and SPT We use the fine-structure transitions of atoms and ions 
as a source of frequency dependent optical depth (r^). These transitions give different contributions to the power spectrum of 
CMB in different observing channels. By comparing results from those channels, it is possible to avoid the limit imposed by 
the cosmic variance and to extract information about the abundance of corresponding species at the redshift of scattering. For 
^■f^ ' PLANCK HFI we will be able to get strong constraints (10"* - 10"^ solar fraction) on the abundances of neutral atoms like C, 

, O, Si, S, and Fe in the redshift range 1-50. Fine-structure transitions of ions like CII, Nil or OIII set similar limits in the very 

important redshift range 3-25 and can be used to probe the ionization history of the universe. Foregrounds and other frequency 
dependent contaminants may set a serious limitation for this method. 
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1. Introduction for all those attempts to estimate, with extremely high preci- 
sion, the key parameters of our universe using this sensitivity 

High precision observations of CMB anisotropies are giving ^^^^^ Efstathiou & Tegmark 1997, Einsestein, Hu & Tegmark 

^ us unique information about the angular distribution of CMB ^^^^^ p^^^^^^^ ^ 2OOO). Indeed, after the first 12 

■ - fluctuations, as well as their spectral dependence in a very f .■ •nT\/tAn . • • »u <; * 

f f J months of operation, the WMAP team is recovering the first 

broad frequency range. HFI and LFI detectors of PLANCK ,.. t t^u r^n/t-o » -^u t 

^ . ? .... multipoles of the CMB power spectrum with an accuracy of a 



spacecraft will provide unprecedented sensitivity in 9 broad ^^^^^^^ (Hinshaw et al. 2003) 
band (Av/v ~ 20 - 30%) channels, uniformly distributed in 
the spectral region of the CMB where contribution of differ- 
ent foregrounds are expected to be at a minimum. CMBpol 



In this paper we are presenting an additional use of the 
tremendous sensitivity of PLANCK and CMBpol, and ground- 
based experiments like ACT, APEX and SPT. We propose to 



and other proposed missions are expected to reach noise lev 

, ^„ .„„ . , ^, ^, ^ „riT AXT^r^iir-T -.u . u look for Or to placc Upper limits On the abundances of heavy el- 

els 20 - 100 times lower than that of PLANCK HFI with tech- t- ft- j 

, 1 J 111 //-.r. 1. '-irvrv'-i\ rr, ^^,^J*T^2 . ii-. cments present in the inter galactic medium and/or in Optically 

nology already available (Church 2002). The WMAP^ satellite ^ ° 1 t- j 



thin clouds of gas everywhere in the redshift range [1,500]. 
We shall focus on the fine-structure lines of neutral (CI, OI, Sil, 



was designed to provide measurements of the CMB temper 

ature anisotropies in the whole sky with an average sensitiv „ , ^ , . . , ^ , . , . , 

j^-ic-T^ rv-io r,io -i i j r , ■ ■ Fel, ...) and lonizcd (CII, NIL OIII, ...) atoms, which might pro 
ity of 35//K per 0.3 x 0.3 pixel at the end of the mission • ^ „ u u , , r ^ a 

(Bennett et al. 2002, Page et al. 2002). Once such sensitivity 

limit is reached, the WMAP data wiU be the first real-Ufe test 



vide information about the epoch of first star formation and ion- 
ization history of the universe. Limits on abundances of heavy 
atoms and ions can be obtained by utilizing the frequency de- 
Send offprint requests to: K.Basu, pendent opacity that will be generated by scattering of back- 

e-mail: kaustuv@mpa-garching.mpg.de ground photons by these species; giving rise to different (JC/-S 

' PLANCK URL site: http://astro.estec.esa.n1/Planck/l in in different observing bands of the experiments. Although 

^ WMAP URL site: http://map.gsfc.nasa.gov/ ^ Planck HFI is expected to detect the combined signal from the 
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strongest lines, like OI 63.2//, CII 157.7// and OIII 88.4/i, fu- 
ture multichannel broadband CMB anisotropy experiments like 
CMBpol might permit to detect contributions from these lines 
separately from the epoch when dark ages were terminating, 
universe became partially ionized and heavy element produc- 
tion begun. It is important that the polarization signal arising 
due to resonance scattering depends strongly on the properties 
of the transition (Sazonov et al. 2002), and together with the 
temperature signal will permit us to separate contributions of 
different species. 

The observed primordial acoustic peaks and angular fluc- 
tuations should not depend on the frequency at all. This is 
connected with the nature of Thomson scattering which pro- 
duces these fluctuations both in the time of recombination of 
hydrogen in the universe (Peebles & Yu 1970, Sunyaev & 
Zel'dovich 1970) and during the secondary ionization of the 
universe (Sunyaev 1977, Ostriker & Vishniac 1986, Vishniac 
1987, and more recently Gruzinov & Hu 1998, da Silva et 
al. 2000, Seljak, BurweU & Pen 2000, Springel, White & 
Hemquist 2001, and Gnedin & JafFe 2001). In this context, 
WMAP polarization measurements have recently shown strong 
evidence favouring an early reionization scenario, with Zr - 
20+^° (Kogut et al. 2003, 95% confidence). But if any amount 
of chemical elements are present during the dark ages, then 
these species will be able to scatter the CMB in their fine- 
structure fines. This scattering would not only partially smooth 
out primordial CMB anisotropics, but in addition will gener- 
ate new fluctuations through the Doppler shift of the line as- 
sociated with the motion of matter connected with the growth 
of density perturbations. The main diff'erence with Thompson 
scattering is that the latter is giving us equal contribution over 
the whole CMB spectrum, whereas the discussed line scat- 
tering would give difi'erent contribution to different observing 
channels placed at different parts of the spectrum. Likewise, 
the contribution in this case would be restricted to a very thin 
slice in the universe. Hence there is a possibility to detect con- 
tribution from the lines with high transition probabilities, even 
though the typical optical depth we would find is very small 
(tv < 10""*). Since every line is able to work only in a given 
range of redshift, having observations on different wavelengths 
can give an upper limit to abundances of different species at 
different epochs. This method relies critically in the fact that 
intrinsic CMB temperature fluctuations are frequency indepen- 
dent. Therefore, in the absence of other frequency-dependent 
components, the difference of two CMB maps obtained at dif- 
ferent frequencies must be sensitive to the difference in abun- 
dance of resonant species at those redshifts probed by the fre- 
quency channels. For this reason, this method is particularly 
sensitive to the possible presence of any frequency-dependent 
signal coming from foregrounds. 

The idea to use the fine transitions of atoms and molecules 
for modifying the CMB power spectrum is not new. Dubrovich 
(1977, 1993) proposed the use of rotational lines of primor- 
dial molecules (LiH, HD etc.) as a source of creating new an- 
gular fluctuations in the CMB, and there have been attempts 
to observe these molecules at high redshifts (de Bernardis 
et al. 1993). Varshalovich, Khersonskii and Sunyaev (1978) 
were trying to estimate the spectral distortions of CMB due 



to absorption in fine-structure lines of oxygen and carbon at 
redshifts around 150-300 where the temperature of electrons 
should be lower than the temperature of CMB due to different 
adiabatic indices of radiation and matter (see Zel'dovich, Kurt 
& Sunyaev 1968). Suginohara et al.(1999) probed the possi- 
bility of detecting excess flux due to emission in these lines 
coming from very over-dense regions in the universe. Recently 
Loeb (2001) and Zaldarriaga & Loeb (2002) (hereafter ZL02) 
computed the distortions connected with the recombination of 
primordial Lithium and scattering in the Lithium resonant dou- 
blet line. Unfortunately the wavelength of fine-structure 2P-2S 
transition of Lithium atoms is too short and will be unobserv- 
able by PLANCK and balloon instruments. On the very low 
frequency domain, future experiments like SKA and LOFAR 
might detect signal from neutral hydrogen 21 cm line, which 
also carry important ioformation from high redshifts (Sunyaev 
& Zel'dovich 1975, Madau, Meiksin & Rees 1997). We con- 
sider below the ground-state fine-structure lines of heavier ele- 
ments, with wavelengths of the order of 50 - 200/i, which will 
in principle be observable with PLANCK and CMBpol if they 
are present in the redshift range [1, 500]. The problem of over- 
coming their extremely small optical depth comprises the main 
idea of this paper. 

We will not discuss in detail the origin and ways of enrich- 
ment of the inter-galactic gas by heavy elements. This certainly 
requires existence of massive stars, supernova explosions, stel- 
lar and galactic winds, and even jets from disks around young 
stars with cold molecular gas (Yu, Billawala & Bally 1999). 
The main goal of this paper is to show that the announced sensi- 
tivity of PLANCK detectors might permit us to set very strong 
upper limits to the time of enrichment of inter-galactic gas 
by heavy elements, the time of reionization, and maybe even 
to detect the heavy elements in the inter-cluster medium. The 
census of baryons in the local universe (Fukugita et al. 1998) 
shows that most of the baryon remains unobserved, and the pro- 
posed method might set way to detect its existence at high red- 
shifts, when it had moderate or low temperature. These miss- 
ing baryons are centainly out of stars, interstellar gas and in- 
tergalactic gas in clusters and groups of galaxies. However, we 
know that such baryons should exist because they have been de- 
tected by WMAP at the last scattering surface at z ^ 1 100, and 
are also necessary to justify the observed abundance of deu- 
terium and ^Li in the early universe. Due to the wide range of 
redshifts under study, and the uncertain degree of mixture and 
dumpiness of our species in the interstellar and inter-galactic 
medium, we shall assume that all elements are smoothly dis- 
tributed in the sky. In other words, in the present paper we shall 
address exclusively the homogeneous low density optically thin 
case, where effects related to strong over-density of gas and 
coUisional excitation are excluded and left as subject of an up- 
coming work. Resonant scattering effects will produce the dis- 
cussed signal even if all the gas in the enriched plasma has over- 
density up to 10^. Under this smooth approximation, we will 
find that the effect of resonant species on the CMB power spec- 
trum wiU be particularly simple, especially in the high multi- 
pole range. Indeed, in the small angular scales, we shall show 
that the induced change in the Q's is given by 6Ci ^ -2 t^jQ, 
where tx, is the optical depth induced by the resonant scatter- 
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ing of species X and C/ is the primordial CMB power spectrum, 
which is currently the main target of most CMB experiments. 

In section 2 we describe the basic approach for computing 
optical depths and corresponding deviations in the CMB 
power spectra due to scattering, and then show how we 
can put constraints on abundances using the sensitivity of 
present and future CMB experiments. Section 3 presents our 
main results for atoms and ions, and in section 4 we discuss 
various possible ionization and enrichment scenarios of the 
universe and the nature of predicted signal for each case. 
We will demonstrate in the Appendix B that optically thin 
over-dense regions with < 10 - 200 cm"-' are contributing 
to the effect of resonant scattering. The effect of foreground 
emissions on our analysis is discussed in section 5, and we 
conclude in section 6. Throughout this paper we have used 
the parameters for standard ACDM cosmology (Spergel et al. 
2003), with Ho = 71 km j"' Mpc-\ Qb = 0.044, £!„ = 
0.27, Qa - 0.73, and tj- = 0.17 for secondary reionization. 
To is the present-day value of mean CMB temperature, given 
by To = 2.726K (Mather et al. 1994). 



2. Basic Approach 

We now discuss the basic method of obtaining the deviations in 
the CMB power spectrum by resonance scattering of atoms and 
ions, and using this deviation to constrain their abundance dur- 
ing the dark ages. The interaction of CMB photons with atoms 
and ions will mostly consist of resonant scattering with either 
an atomic or a rotational/vibrational transition, depending on 
the species under study. In the Appendix A we detail how we 
introduce the optical depth due to the resonant transition (tx,) 
in the CMBFAST code (Seljak & Zaldarriaga, 1996). We con- 
sider a particular resonant transition i for a given species X, 
with rest-frame resonant frequency v,. The total optical depth 
encountered by CMB photons on their way from the last scat- 
tering surface to us is then obtained by adding the contributions 
from all Unes to the standard Thomson opacity: 



T = Tt + 



(1) 



In order to calculate tx, , we shall recur to the formula which 
gives the optical depth of a resonant transition in an expanding 
medium (Sobolev 1946), 



TXiiz) = fi 



7Te^ AiUxiiz) 
nieC H(z) 



(2) 



where is the absorption oscillator strength of the resonant 
transition, /I, is the corresponding wavelength (in rest frame), 
nxi(z) is number density of X species at redshift z, and H{z) 
is the Hubble parameter at that epoch. The oscillator strength 
depends on At, the Einstein coefficient of the transition Aj,;, and 
the degeneracy of the levels involved in it: 



fi 



(3) 



A simple and elegant treatment following Gunn & Peterson 
(1965) using 5-function hne profile, instead of thermally broad- 
ened gaussian for the same line, gives same value of optical 



depth. Further modifications are made to this formula to take 
into account the finite populations in the upper transition levels, 
which is important for atomic fine-structure transitions, where 
excitation temperature for them is comparable to the radiation 
temperature at high redshifts. Also the effect of non-zero cos- 
mological constant can not be neglected in the low redshift uni- 
verse of interest. Considering these facts, we obtain the follow- 
ing formula for optical depths: 



rxiiz) 



1.7 X 10" 



10-6 



Ai 



Siz = 10)/ \lOOn) \10 



(4) 



Here Xx,(z) is the ratio between the number density Nxi(z) 
of the atomic or ionic species X under consideration, with 
respect to the baryon number density at the same redshift: 
Nh(z) = 2 X 10-^ (1 + zf cm-l i.e. Xx,(z) = NxJNi, presents 
the evolution of abundance of the given species due to ele- 
ment production, ionization and recombination processes. We 
propose to constrain the minimum abundance that can be de- 
tected at that redshift in solar units, or [X]mi„ = XxXz) I Xq, 
where, for example, = 3.7 x 10"* for carbon. S{z) gives 
the redshift dependence of optical depth in a ACDM universe: 
S{z) = (1 + z? [(1 + zf{l + a„z) - z{2 + z)Qa]4 (see, e.g., 
Bergstrom 1998) , and A, is the wavelength in micron. The fi- 
nal term S accounts for the actual fraction of atom/ion present 
in their ground state, and is governed solely by the tempera- 
ture of background radiation at the redshift of scattering. For a 
two-level system, this fraction is simply 



ni = [l+ igjgi) expi-hvJkToil + zWK 



(5) 



We also include here the correcting term for the induced emis- 
sion in the presence of the CMB and finite population of the 
upper level: 



!B = ni X 



1 - exp - 



hvi 



kBToil + z) 



(6) 



The resonance scattering on ions and atoms in thermal 
equilibrium with black body radiation does not change its in- 
tensity. However, the observed CMB also has finite primordial 
angular fluctuations. The effect of resonant scattering is to de- 
crease these angular fluctuations, and to bring the system more 
close to thermodynamic equilibrium. Therefore resulting fluc- 
tuations observed on the frequency of the line should differ 
from the situation on other frequencies which are far from the 
resonance. 

At low multipoles peculiar motions arising due to the 
growth of large scale density perturbations become important. 
All ions or atoms are moving in the same direction and change 
the frequency of CMB photons during resonant scattering. This 
leads to generation of new anisotropics of background. 

We next give a simple description of the modification of 
the power spectrum of the CMB when it encounters a resonant 
species. As mentioned in the Introduction, our first hypothesis 
will be that the species responsible for the scattering are homo- 
geneous, isotropic and smoothly (i.e., not clumped) distributed 
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in the Universe, at least during the epoch of interaction with the 
CMB. One may argue that heavy species are located in halos 
where the first stars are produced and that their distribution in 
the sky can non be regarded as smooth. In this case, the final 
effect would depend on the typical angular size of the patches 
in which the species have been spread and on their total sky 
coverage. This paper will observe the case where those scatter- 
ing patches percolate in the sky, giving rise to a smooth, ho- 
mogeneous picture. The case of patchy distribution of emitting 
sources will be addressed in a forthcoming work, where col- 
lissional processes are studied in an extremely dense optically 
thick environment. 

In the conformal Newtonian gauge, (also known as the lon- 
gitudinal gauge), the A;-mode of the photon temperature fluctu- 
ation at current epoch is given by (Hu & Sugiyama 1994) 

Xm 
dr^e'''^'^''-'^'>\r{r^l){^To-i^iVb) 

+ ^ + i// - ikfiiffj. (7) 



Here we have neglected the polarization term (which con- 
tributes at most with a few percent of the temperature ampli- 
tude). T] denotes the conformal time, Y = t(?7) exp -t(?7) is the 
visibility function and t(?7) - J^" drj'firj') is the optical depth. 
Aro accounts for the intrinsic temperature fluctuations, vi, for 
the velocity of baryons and (p and i/r are the scalar perturbations 
of the metric in this gauge. We have also neglected all tensor 
perturbations. If we now introduce the optical depth associated 
to a resonant transition as a Dirac delta of amplitude rx, placed 
at// = Tlx,, 



dri?]) = o-TrieirOdri + Tx,6(ri - rix,)dri, 



(8) 



we readily obtain that the original anisotropics have been 
erased by a factor e '^^f , whereas new anisotropics have been 
generated at the same place: 



A7-(fc,?7o,yu) = e '"'iAT„,Jk,rix,,fj) + AT„Jk,T]x,,fj) 



(9) 



In real space, if we assume that <k 1 , the Dirac delta ap- 
proximation for the resonant transition (eq.(|S|l) translates into 

— (n,z = 0) = (l-Tx,) — (n,zx,) 

JO ^0 



5T 



lin 

(n,zx,) + 0[tI]. 



(10) 



Here, n is the unitary vector denoting the observing direction, 
and 5r/ro|J"„,(n, zx,) is flie linear in t^, term of the newly gen- 
erated anisotropics-'. We remark again that we are assuming 
that the species is homogeneously distributed, so that txi does 
not depend on position. Let us now consider two dififerent ob- 
serving frequencies: the first one corresponds to the scattering 
redshift zx,, and the another corresponds to a scattering redshift 
too high to expect any significant presence of the species X. If 
we compute the correlation function, ^||-(ni, z = 0) ■ f|-(n2, z = 

^ In what follows, our temporal coordinate will be denoted by red- 
shift or conformal time r] indiferently. 



0)^, in both cases, we find that the difference of these quantities 
will be equal to:"* 



-2(— (ni,zx,)- — (n2,zx,)) 



+ 2( — (ni,zx,)- — 



(n2,zx,)))+(5[T|J. (11) 



That is, the term linear in tx, of 
A((ff(ni,z = 0)- fS(n2,z = 0))) is the sum of two dif- 
ferent contributions: the suppression of original fluctuations 
(which, as we shall see, dominates at small angular scales), and 
the cross-correlation of the newly generated anisotropics with 
the intrinsic CMB field. The first term (suppression) does not 
depend on the potential or velocity fields during scattering, but 
only on the intrinsic CMB anisotropy field. Moreover, we must 
remark that those contributions are evaluated at scattering 
epoch^. Hence, if the resonant scattering takes place before 
reionization, the changes in the correlation function (or in the 
power spectrum, as we show below) will be sensitive to the 
CMB anisotropy field before being processed by the re-ionized 
medium. In the Appendix A we give a detailed computation 
of the change in the power spectrum due to a resonant line. As 
we concentrate in the optically thin limit, it is possible to make 
a power expansion on the optical depth tx, : 



5Ci = Cf - C, = Tx, ■ Ci + 4 ■ Ci + 0(tI). 



(12) 



In this equation, C, ' and C; refer to the angular power spec- 
trum multipoles in the presence and in the absence of the X, 
resonant transition, respectively. For the limit of very small 
Tx, one can retain only the Unear order, and a direct identifi- 
cation of 5Ci with the abundance {X\{z) is possible by means 
of eq.©: the change induced in the CMB power spectrum will 
be proportional to the abundance of the species responsible for 
the resonant scattering. In agreement with what has been es- 
tabUshed when studying A (^|?(iii,z = 0) ■ f^(n2,z = 0)^), we 
find that C I changes sign at some multipole Icwss'- positive val- 
ues of 6Ci imply generation of new anisotropics {I < Icross), 
whereas anisotropics are suppressed for I > Ic-oss (see figure ^ 
at very low I). In the Appendix A we note that positive values 
of 6Ci (generation) are due to the non-zero correlation exist- 
ing between anisotropics generated during recombination and 
those generated during the resonant scattering. This correlation 
is due to the coupling of k modes of the initial metric pertur- 
bation field in scales of the order or larger than the distance 
separating the two events (recombination at redshift 1100 and 
resonant scattering at redshift 5-25), (Hernandez-Monteagudo 
& Sunyaev, in preparation). 

This type of manipulations of the correlation function/power spec- 
trum which avoids the Cosmic Variance limit is possible as long as the 
weak signal (signal induced by scattering in this case) has different 
spectral dependence than the CMB. 

^ The dependence on cosmic epoch of eqs. 11011 it has been simpli- 
fied. We again refer to Appendix A for a formal derivation. 
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Fig. 1. The nature of temperature anisotropy that will be 
marginally detectable with Planck HFI 143 GHz channel, using 
63 fi OI line. Abundance is taken from the tableQ, as 5.3 x 10""^ 
solar at redshift 32. Here Primary denotes the measured tem- 
perature anisotropy (upper line), C'l"^' + 6Ci ~ C'l" , and Line 
denotes the the newly generated anistropy, |i5C/|, arising from 
line scattering {bottom, with filled squares). This \6Ci\ is ob- 
tained by taking the difference from the 100 GHz channel, and 
changes sign around Z = 3. The corresponding noise level is 
denoted by HFI noise. Also shown is the cosmic variance limit 
{short-dashed line) for comparison. 



2.1. dCi's at small angular scales 

With respect to the high-1 range, we find that the change in- 
duced in the power spectrum takes a very simple form: 



6C, - -2 Tx. Ci 



(13) 



where Q is the intrinsic power spectrum. This dependence 
is identical to the effect of reionization on the power spectrum 
at small angular scales. Indeed, in that scenario, if the optical 
depth due to electron scattering during this epoch is given by 
T,.f,o, then we have that for / » 1 the intrinsic CMB power 
spectrum generated at recombination is suppressed by a fac- 



tor 1 - e.yip{-2Treio), or 



-2t/-^(V) if 



< 1. So in both 



(resonant and Thompson) scatterings, the shape of the induced 
change in the power spectrum is particularly simple and equal 
to 6Ci ^ -2t Ci, for a given optical depth t and intrinsic power 
spectrum C;. This simplifies considerably the effect of reion- 
ization on the dCi's induced by resonant transitions. Indeed, if 
the symbol A'''^'" denotes the difference of a given quantity eval- 
uated in the presence and in the absence of reionization, then 
we have that, for high /, A"''''(C/) = -iTreu, Ci and A'"" {6Ci) = 
-iTreio {dCi). This is shown in figure 0, where we plot the 
quantities V = \Cf°ICi - l| and W = \5Cf°l5C\ - l|, that is, 
the relative change of C\ and 6C\ due to the presence of reion- 
ization. We have taken t„,o = 0.17. We see that, for high /, both 
y and W approach the limit y ^ W ^ 2 Treio^ i-c, the effect of 
reionization is identical in the two cases. 




Fig. 2. For two lines placed at different redshifts, (z = 20 in 
the top, z - 10 in the bottom), we plot the quantities Y (solid 
lines) and W (dotted lines), defined as y = \C\"" ICi - l| and 
W = ||(5C['"°/(yC;| - l|, respectively. In both cases, in the limit 
of high these quantities tend to 2 T,-f,o, 
to be 0.17. In both panels, tx, - 10"^ 



2.2. Measuring 6Ci 's and abundances 

Next we investigate the limits of detectability of the optical 
depth of atoms and ionss in current and future CMB missions. 
Our starting point is the expected uncertainty in the obtained 
C;'s from any CMB experiment (Knox 1995): 



(2/+1)/, 



skv 



(14) 



In this equation, C; is the underlying CMB power spectrum, 
fsiiy is the fraction of the sky covered by the experiment, and 
w is the pixel weight given by w"' = cr^fi;„v, with cr^ the 
noise amplitude for pixels of solid angle Q.pix- Bi is the beam 
window function, and we shall approximate it by a gaussian. 
We are assuming that the noise is gaussian, uniform and un- 
corrected. The first term in parentheses of ea.( ll4t reflects the 
uncertainty associated to the cosmic variance. In the ideal case 
of an experiment with no noise, this would be the unavoidable 
uncertainty when identifying the observed C/'s to the C;'s cor- 
responding to a particular cosmological model. However, our 
interest focuses on the comparison of power spectra measured 
in different frequency channels with respect to a measured ref- 
erence power spectrum, (which is supposed to be free of any 
contamining species of atoms and/or ions). Under the assump- 
tion that the reference power spectrum contains only CMB and 
noise components, we can write the following expression for 
the uncertainty in the measured power spectrum difference be- 
tween a probe channel and the reference channel: 



'SCi 



{2l+\)fsky 



prob 



-1 

prob^ l,prob 



.B 



(15) 



The indexes prob and re/ refer to the probe and reference chan- 
nel respectively, and the first term dC^"''^ refers to the cosmic 
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variance associated to the temperature anisotropics generated 
by the resonant species. 

Let us now assume that a CMB experiment is observing 
at a frequency for which it is expected to see the effects of a 
resonant transition of a species X at redshift 1 + zx, = vxjvots- 
One can then make use of ea. J12> to relate the minimum tx, 
observable with the sensitivity of the experiment: 

The factor n expresses the cr-level necessary to claim a detec- 
tion. For all the results presented in this paper, we have set the 
detection threshold at 3cr, i.e. n = 3. Fig[n shows how a very 
low optical depth suffices to generate changes in the power 
spectrum (filled squares) which are big enough to overcome 
the combined noise of the two experiment channels being 
compared, (dotted line). Once this minimum optical depth 
is obtained, we can easily find the corresponding minimum 
abundance through eq.Q, which is one of the main goals 
of our paper. In fact, this limit can be further improved by a 
factor of ~ VAZ if one computes band power spectra on some 
multipole range Al. We shall address the details of these issues 
in the next section, where we show that for PLANCK' s HFI 
channels, values for the optical depth as low as 10"^ - 10"^ 
can be detected. 

So far we are neglecting the problematic associated to 
the calibration between channels and the possible presence 
of foregrounds. Given the frequency dependent nature of 
the latter, the amplitudes and their characterization of the 
different galactic and extragalactic contaminants are critical 
for our purposes. We return to these issues in section 5, were 
we estimate the eff'ect of foreground contamination on our 
analysis. 



3. Results and Discussion 

Now we present the results for selected atoms and ions and that 
can produce measurable distortions in the CMB spectrum. As 
mentioned previously, the focus will be on the fine-structure 
transitions for atoms and ions , as their frequencies fall in the 
far-infrared and microwave range and are therefore perfectly 
suited to our purpose. 

We divide this section into two parts: the first part shows 
results for all the important atomic and ionic species. In sec- 
ond part we briefly mention the contribution from the ovedense 
regions to the effect considered. 

3.1. Scattering by atoms and ions of heavy elements 

In this section we investigate the possibility of distorting the 
CMB power spectrum by scattering from neutral atoms like CI, 
OI, Sil, SI, Fel etc., as well as singly and doubly ionized species 
of heavy elements, like CII, Nil, Sill, Fell, OIII etc. This is im- 
portant for various reasons. The Gunn-Peterson effect permits 
us to prove that the universe was completely ionized as early 
as redshift 6, up to the position of most distant quasars known 



today. Recent results from WMAP satellite pushes the reioniza- 
tion redshift as far as z = 20, suggesting a complex ionization 
history (Kogut et al, 2003). There is extensive discussion about 
the nature of reionization, and also on the possibility for uni- 
verse being reionized twice (Cen 2003). In any case, before the 
universe was ionized completely, there were regions of ionized 
medium around first bright stars and quasars. It will be diffi- 
cult to prove that the universe was partially ionized using Ly-a 
line because of its extremely high oscillator strength, which 
makes the universe optically thick in this line even when the 
neutral fraction of Hydrogen is only 10 at z = 6. However, 
the infrared lines that we are discussing in this paper have much 
weaker oscillator strength and therefore neutral gas is transpar- 
ent in these lines up to very high redshifts, even if we assume 
solar abundance. This fact permits us to consider the possibil- 
ity that Planck and other future CMB experiments setting very 
strong limits on abundance of neutral atoms in the early uni- 
verse. 

The effect of line emission to the thermal spectrum of CMB 
can be estimated in a simple order-of-magnitude way by the 
formulation given in Appendix B. If there were significantly 
over-dense regions at high redshifts (up to z ~ 20) which were 
completely ionized and enriched with metal ions (e.g. C^, N^, 
0++, Fe^ etc.), then collisional excitation followed by radia- 
tive de-excitation will be a significant source of emission in the 
same fine-structure lines. However, to make line emission visi- 
ble we need three factors: high abundance of the elements, high 
density of the electrons in the strongly over-dense regions and 
large amount of over-dense regions in the volume which we are 
investigating. Hence it is much more promising to study angu- 
lar distortions of the CMB generated by scattering from the low 
density regions of the universe, rather than studying distortions 
in its thermal spectrum, for constraining heavy element abun- 
dances at high redshifts. Such low density inter-galactic gas is 
believed to contain most of the baryonic mass of the universe, 
and possibly exists as warm/hot gas with 10^ < T,. < 10^ K 
today (Cen & Ostriker 1999). But at redshifts z > 1 this inter- 
galactic gas should have moderate (T^ ~ 10"* K) or low (T^, ~ 
Tcmb(z)) temperature, and the proposed method of observing 
angular fluctuations caused by scattering from neutral or singly 
ionized atoms might set a direct way to detect its existence. 

As mentioned above, the main contributors of opacities 
in the relevant frequency range are oxygen, nitrogen, carbon, 
sulfur, silicon, and iron, along with minor contributions mainly 
from phosphorus, aluminum, chlorine and nickel. The 63// 
fine-structure line of neutral oxygen gives strong constraint on 
neutral species at high redshift, but early reionization makes 
lines of CII, Nil and OIII even more important. All data 
relating to fine-structure lines have been taken from the ISO 
line-list for far-IR spectroscopy (Lutz, 1998), and ffie Atomic 
Data for the Analysis of Emission Lines by Pradhan and Peng 
(1995). When necessary, this compilation was supplemented 
by the freely available NIST Atomic Database*". 

The basic idea of obtaining limits on abundances can be 
understood from Fig^and tabled Each broad-band channel 

^ http://physics.nist.gov/cgi-bin/AtData/Iines_form 
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Atom/ 


Wavelength 


Oscillator 


HFI freq. 


Scattering 


S 


Opt. depth for 


[X]„„„ for 


<M,m„> in 


Ion 


(inyu) 


strength 


(GHz) 


redshift 


factor 


IQ-^ solar abundance 


/ = 10 


/ = 10 - 20 


CI 


609.70 


1.33 X 10-" 


143 


2.4 


0.76 


6.4 X 10-^ 


5.3 X 10-' 


2.6 X 10-" 








217 


1.3 


0.92 


3.9 X 10-"^ 


1.4 X 10-- 


6.8 X 10-" 








353 


0.4 


0.99 


1.6 X 10-^ 


2.1 X 10-' 


1.2 X 10-' 




370.37 


9.08 X 10-'" 


143 


4.7 


0.15 


1.2 X 10-^ 


2.8 X 10-- 


1.3 X 10-- 








217 


2.8 


0.09 


3.7 X 10-^ 


1.6 X 10-' 


8.1 X 10-^ 


CII 


157.74 


1.71 X 10-" 


143 


12.3 


0.79 


1.8 X 10-^ 


2.7 X 10-- 


6.2 X 10-" 








217 


7.9 


0.94 


1.1 X 10-5 


7.7 X 10-^^ 


3.0 X 10-" 








353 


4.4 


0.99 


5.6 X 10-*^ 


7.7 X 10-- 


3.6 X 10-- 


Nil 


205.30 


3.92 X 10-" 


143 


9.2 


0.76 


1.1 X 10-5 


7.6 X lO-"* 


2.6 X 10-" 








217 


5.8 


0.92 


6.8 X 10-"^ 


8.6 X 10-^^ 


3.8 X 10-" 








353 


3.1 


0.99 


3.5 X 10-'^ 


1.3 X 10-' 


6.8 X 10-- 




121.80 


2.74 X 10-" 


143 


16.2 


0.16 


2.1 X 10-^ 


1.3 X 10-' 


3.8 X 10-- 








217 


10.5 


0.09 


6.4 X 10-'' 


3.4 X 10-' 


1.1 X 10-' 


NIII 


57.32 


4.72 X 10-" 


143 


35.6 


0.79 


2.5 X 10-5 


2.3 X 10-3 


7.4 X 10-* 








217 


23.4 


0.94 


1.5 X 10-5 


6.1 X 10-^^ 


2.0 X 10-" 


OI 


63.18 


3.20 X 10-" 


143 


32.2 


0.88 


1.0 X 10-* 


5.3 X 10-* 


1.7 X 10-* 








217 


21.2 


0.96 


6.3 X 10-5 


2.0 X 10-' 


6.4 X 10-* 








353 


12.5 


1.00 


3.1 X 10-5 


2.2 X 10-' 


4.9 X 10-- 


OIII 


88.36 


9.16 X 10-" 


143 


22.8 


0.76 


2.2 X 10-* 


3.5 X 10-* 


1.2 X 10-* 








217 


14.8 


0.92 


1.4 X 10-* 


8.4 X 10-3 


1.8 X 10-" 








353 


8.6 


0.99 


7.4 X 10-5 


1.2 X 10-- 


4.4 X 10-" 




51.81 


6.55 X 10-" 


143 


39.5 


0.17 


4.5 X 10-5 


1.4 X 10-" 


4.7 X IQ-* 








217 


26.0 


0.10 


1.4 X 10-5 


6.5 X 10-3 


2.2 X 10-" 


Si I 


129.68 


6.24 X 10-" 


143 


15.2 


0.76 


6.4 X 10-^ 


6.9 X 10-- 


1.8 X 10-- 








217 


9.8 


0.92 


4.2 X 10-* 


3.7 X 10-- 


1.2 X 10-- 




68.47 


4.92 X 10-" 


143 


29.7 


0.19 


1.8 X 10-* 


3.1 X 10-- 


1.0 X 10-- 


Sill 


34.81 


7.74 X 10-" 


217 


39.2 


0.94 


2.0 X 10-5 


4.8 X 10-" 


1.6 X 10-" 








353 


23.4 


0.99 


1.0 X 10-5 


7.4 X 10-" 


2.4 X 10-- 


SI 


25.25 


8.03 X 10-" 


217 


54.4 


0.96 


6.0 X 10-* 


2.4 X 10-- 


7.4 X 10-" 


Fel 


24.04 


1.69 X 10->* 


143 


86.4 


0.83 


3.2 X 10-5 


4.9 X 10-" 


1.3 X 10-" 








217 


57.2 


0.95 


2.0 X 10-5 


7.7 X 10-" 


2.3 X 10-" 








353 


34.4 


0.99 


9.7 X 10-* 


7.5 X 10-- 


2.5 X 10-- 




34.71 


2.06 X 10-"* 


143 


59.5 


0.02 


9.2 X 10-' 


1.1 X 10-' 


3.4 X 10-" 


Fell 


25.99 


1.73 X 10-'* 


217 


52.9 


0.95 


1.9 X 10-5 


7.1 X 10-" 


2.2 X 10-" 








353 


31.7 


0.99 


1.3 X 10-'' 


7.4 X 10-- 


2.4 X 10-- 


Felll 


22.93 


3.09 X 10-'* 


353 


36.1 


0.99 


1.8 X 10-5 


4.2 X 10-- 


1.4 X 10-- 



Table 1. Minimum abundance of the most important atoms and ions that can be detected from Planck HFI. We have used a fixed 
value of abundance (1% solar) to obtain the optical depth in the fine- structure transitions of any given species, and because of the 
very low values of optical depths, have used the linear relation from eqn.(ll6ll to obtain the minimum detectable abundance for 
the sensitivity limit of Planck HFI. The S factor is the correction term obtained from eqn.(|6j. In the last column we have further 
improved the sensitivity by averaging the instrument noise in the multipole range / = 10 - 20. The 100 GHz channel has been 
used as reference for all the cases. We present central redshift for coiTesponding channel, however, in reality HFI will be able to 
give limits only for redshift intervals corresponding to the widths of frequency channels. 



of HFI acts in a specific range of redshifts for a particular line, 
and we have tabulated the central redshift corresponding to the 
scattering for three most sensitive channels for several atomic 
and ionic fine-structure transitions. The lowest frequency chan- 
nel of 100 GHz is assumed to be "clean" from line scattering, 
and thus used as reference. We use a fixed abundance (lO"- 
solar) to obtain the optical depth in accordance with formula 
0. Such small values of optical depths allow us to use the 
first-order approximation, and so we finally obtain the mini- 
mum optical depth , and hence the minimum abundance (last 
column, with respect to solar value) from the sensitivity level 
of the detector (at 3cr level). We have neglected signals below 
/ = 5, especially at the quadrupole or Z = 2 where noise level is 



minimum, due to the fact that it will be very difficult to observe 
the predicted signal at such large angular scales due to the fore- 
grounds. The best angular range for Planck HFI is / ~ 10 - 30, 
and we can average it over a multipole range Al and improve 
the detectibility by a factor of VA/. This result is shown in the 
last column of our table. 

Fig^ shows the expected behavior of (5C/-S when neutral 
oxygen is marginally detectable with the 143 GHz channel of 
HFI, limiting the oxygen abundance as low as 10"^ relative to 
solar at redshift 30. We show both the measured temperature 
anisotropy and 6Ci-s generated by line scattering in this plot. 
This line contribution touches the noise limit at around I ~ 10, 
showing the angular scale where best possible limit can be ob- 
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tained. We see that this effect always Hes much below the cos- 
mic variance limit, but due to the frequency dependence of new 
anisotropies we are not constrained by this limit. 

Part of the results of this paper are summarized in Tabled 
This lists all the atoms and ions on whose abundances we can 
put strong limits, and all these limits are computed at 3cr level. 
With HFI we have the possibility to use more than one probe 
channel to show different upper limits for the same species at 
different redshifts. This fact can be helpful to model the abun- 
dance history of the universe, and we present a general discus- 
sion in the next section. 

3.2. Contribution from over-dense regions 

Enrichment of primordial gas by heavy elements occurs due to 
supernova explosions of the first stars. High velocity stellar and 
galactic winds and low velocity jets from disks around forming 
stars, and objects of the type of SS 433 with baryonic jets carry 
enriched matter to a large distance from the forming stars. The 
observation of the most distant galaxies and quasars are show- 
ing that even most distant objects (z ~ 5 - 6) have chemical 
abundance on the level of solar (Freudling at al. 2003, Dietrich 
et al. 2003). At the same time there is a possibility that the 
low density matter, e.g. in future voids, will have extremely 
low abundance of heavy elements. Observational method we 
are proposing might permit us to observe ions and atoms of 
heavy elements in diffuse matter, with over-density lower than 
10^- lO'* at redshifts 10-20, and up to 5 ~ 10^- 10*' at z ~ 2-5. 
This means that even Ly-or clouds are contributing to our effect. 
We should be careful only with the most damped Ly-cc systems, 
and with dense gaseous nebulae of the type of Orion and dense 
giant molecular clouds. Diffuse gas in the galaxies and proto- 
galaxies should also contribute to our signal. In Appendix B we 
estimated the level of densities in the gas clouds, when the dis- 
cussed effect is diminished by 30% or 50% in comparison with 
the case of diffuse inter-galactic space, and showed that over- 
densities greater than 10^ are needed at high redshifts before 
collisional effect begins to decrease amplitude of scattering sig- 
nal. To summarize, the proposed method might permit us to get 
signal from all diffuse matter of the universe, excluding only 
the extremely over-dense, or optically thick clouds. It is proba- 
ble that the over-dense regions were the first to be enriched by 
heavy elements. Therefore it is important that the moderately 
over-dense regions of the universe are also contributing to the 
resonance scattering signal in the CMB angular fluctuations. 

The simultaneous effect of free-free, line and dust emis- 
sion from non-uniformly distributed over-dense regions at high 
redshifts, which are in non-linear stage of evolution and enter- 
ing the state of intense star formation, produces an independent 
signal from the scattering effect considered in this paper. The 
same lines which contribute to the C/-s due to resonant scatter- 
ing from extended low density regions, would also contribute to 
the power spectrum at smaller angular scales due to emission 
from over-dense regions. But over-dense objects like damped 
Ly-Q- absorbing systems and low density Ly-c clouds, together 
with stars, atomic and molecular gas in galaxies and hot gas in 
clusters, contains only 20% - 40% of the baryons in the uni- 



verse (Fukugita et al. 1998, Penton et al. 2000, Valageas et al. 
2001); meaning that resonance scattering from low density, op- 
tically thin gas with low temperature and moderate (<10^-10^, 
see Appendix B) over-density will always create its own distor- 
tion in the CMB power spectrum alongside the emission gener- 
ated from denser regions. Both effects carry information about 
the abundances of atoms and ions averaged over the volume 
defined by multipole / (or angular scale), and the frequency 
resolution of the detector which gives us the thickness of the 
slice in redshift space along line of sight. Resonant scattering 
effect is sensitive only to the mean density of the scattering 
species (n^, ) in that volume, whereas the line radiation effect is 
connected with collisions and therefore its contribution to C/'s 
depends on in the same volume. They are, thus, two in- 

dependent effects that carry complementary information. This 
latter effect is sensitive to the most over-dense regions in the 
universe, and will be discussed in detail in a subsequent work 
(Basu, Hernandez-Monteagudo and Sunyaev, in preparation). 

4. The Ionization i-iistory of The Universe 

Detection of the OI, CII and OIII lines of the two most abun- 
dant elements will permit not only to trace the enrichment of 
the universe by heavy metals, but also might open the way to 
follow the ionization history of the universe. According to re- 
cent models of stellar evolution very massive Pop III stars effi- 
ciently produce heavy metals like oxygen, carbon, silicon and 
sulfur (Heger & Woosley, 2002). The CNO burning phase of 
the stars appearing immediately afterwards will also produce 
large amount of nitrogen, and we can expect strong signal from 
the ionized nitrogen 205;U line, as the time of evolution of the 
first stars is extremely short (~ 10^ years) in comparison with 
the Hubble time even at redshift 25. 

The WMAP finding that universe has rather high optical 
depth due to secondary ionization at redshift Zr - 20^^" (Kogut 
et al. 2003, 95% confidence) forced many theoretical groups 
to return to the picture of early ionization due to Pop III stars 
(Cen 2003, Wyithe & Loeb 2003). One of the possible evolu- 
tion scenario of abundances for the elements produced by these 
massive hot Pop III stars and intense star and galaxy formation 
is given in Fig. 3a. Here we consider two enrichment histories of 
the universe, with low (A) and high {E) metal abundances after 
reionization. The first phase of metal enrichment occurs during 
the epoch of Pop III stars, and at later epochs (redshift 3-5) 
intense galaxy formation causes further rise in the metallicity. 
Later in this section we consider a third enrichment history with 
late reionization and metal injection. 

In Fig. 3b and 3c , we sketch the relative ion fraction of the 
three most important atomic and ionic species under above- 
mentioned ionization history. We show two different reioniza- 
tion scenario: one for relatively cold stars when production OIII 
is less efficient (case 1), and the other for hot stars and quasars 
which are able to keep oxygen fully ionized at all intermedi- 
ate redshifts (case If). Line of OI should give us an information 
about creation of oxygen before Universe was strongly ionised. 
Relative growth of CII line (see Fig. 2b) will mark the time 
when carbon will be ionised in large ionised regions which do 
not completely overlap or will be partially ionised everywhere 
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Fig. 5. Angular dependence of temperature anisotropy gener- 
ated by scattering from CII 158// fine- structure line. Shown 
here are three cases of \6Ci\ as should be observed from Planck 
HFI 143 GHz (scattering at z = 12), 217 GHz (scattering at 
z - S) and 353 GHz (scattering at z = 4) channels, (using 100 
GHz channel as reference)for 10% abundance of CII ions (with 
respect to solar) at all three redshifts. Also shown are the pri- 
mary anisotropy, C/, the cosmic variance limit, and noise level 
for 217 GHz channel (with respect to 100 GHz) for compari- 
son. 



(we can not dinstinguish this two variants of ionization his- 
tory using large angle observations). Potential of CI ionization 
(1=11.26 eV) is lower than that of HI. Therefore CII fraction 
might be higher than that of HII and Oil in the beginning of sec- 
ondary ionization. At the same time OI fraction should follow 
that of HI because ionization potentials are so close (1=13.62 
and 13.60 eV correspondently). Pop III stars should be very 
hot, and therefore they are able to ionize helium early enough. 
Simultaneously OIII should become abundant ion because ion- 
ization potential of Oil (1= 35. 12 eV) is higher than that of Hel 
but smaller then that of Hell. 

Figl^shows the amplitude of the predicted signal as a func- 
tion of observing frequency, for a fixed angular scale I - 10. 
The contributions from four most important lines, viz. CII 
158yu, Nil 205yL(, OI 63yL( and OIII 88//, are shown, along with 
their sum. We also show four diff'erent redshift ranges for each 
line, which highUghts the fact that contributions from CII and 
Nil are higher because their signal is coming from lower red- 
shifts where abundance is higher according to our chosen abun- 
dance history. In Fig|5]we show the angular dependence of the 
temperature anisotropy generated by resonant line scattering 
for the three important HFI channels. The abundance of CII ion 
is kept fixed at 10% solar, so that the lower observing frequency 
(i.e. higher scattering redshift) has a higher value of optical 
depth in the same line. The left-side of each curve is dominated 
by Doppler generation of new anisotropy, and hence are posi- 
tive. The righ-side is dominated by suppression of primordial 
anisotropy and hence SCi-s are negative, the discontinuity in 
each curve shows the interval where 6Ci changes sign. We see 
at low multipoles both generation and suppression term tend 
to cancel each other. Especially due to the adoption of WMAP 



reionization model with high optical depth (zreio - 0.17) in 
our computations, lines which scatter CMB photons at red- 
shifts z < 18 encounter a high value of visibility function due 
to reionization, and therefore causes a strong Doppler genera- 
tion (but negative) in the multipole range Z = 20 - 100. This is 
the cause for higher amplitude of the signal at 143 GHz. But 
at high multipoles Z > 100 the contribution from Doppler gen- 
eration is negligible, and the signal is simply proportional to 
the primordial C/-s, and the line optical depth, provided t «: 1 . 
Hence using data from tableQ, and knowing the amplitude of 
primordial CMB signal, one can immediately predict the am- 
plitude of the effect at small angular scales using formulalfT^. 

Fig|6l summarizes our results for three different angular 
scales. As explained above, at low multipoles the effect is not 
proportional to C/, as can be seen from the two panels at top. 
These low multipole results are important for future satellite 
missions like Planck and CMBPol. In the middle, results for 
/ = 220 are shown, where effect is already proportional to C;. 
The amplitude of the signal is highest here due to the large 
amplitude for the C; at the first Doppler peak. At the bottom 
we give results for the third Doppler peak, or Z = 810, where 
signal drops by a factor of 2 from that of / = 220. However, 
these angular scales {0 ~ 13') are particularly suitable for fu- 
ture balloon and ground-based experiments with multi-channel 
narrow-band recievers. 

Finally, in Fig0 we present an alternative ionization and 
enrichment history of the universe, when there were no produc- 
tion of heavy elements before z ~ 12. This is very interesting 
because, even for such late reionization, if there is moderate en- 
richment of the IGM around redshift 4-5, we have the possibil- 
ity to detect our signal with Planck HFI around / = 10. The con- 
tributions from high energy oxygen lines are reduced because 
of the absence of metals at high redshift, but long-wavelength 
lines of CII and Nil still can generate strong signal from low 
redshifts (z < 8). The HFI sensitivities in this figure have again 
been improved by averaging the noise in a multipole range of 
AZ = 10 around / = 10. As discussed before, the signal for 
/ = 220 and / = 810 are proportional to the optical depths in 
these lines and the primordial C/-s. 

Due to relatively small scattering cross-sections of the fine- 
structure transitions under discussion, such observations are 
sensitive to significant abundances of the atoms and ions (for 
example when given species contributes from 10% up to 100% 
of the corresponding element abundance), whereas the Gunn- 
Peterson effect gives optical depth of the order of unity al- 
ready when abundance of neutral hydrogen is of the order of 
10"^ - 10"^. Hence even if the universe is completely opaque 
to Ly-a photons at redshifts z > 6, ionic fine-structure lines 
like CII 157 J fi can probe the very first stages of patchy reion- 
ization process. Detection of all three broad spectral features 
connected with strongly redshifted OI, CII and OIII microwave 
lines will open the possibility to trace the complete reionization 
history of the Universe. In addition OIII line will proof the ex- 
istence of extremely hot stars at that epoch. 

The published level of noise of the Planck HFI shows that 
the three low frequency channels of HFI are almost an order of 
magnitude more sensitive than signal coming from cases I & II 
of history B. Even in the case of late reionization history (case 
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C, Fig0 ), Planck is about 4-times more sensitive than pre- 
dicted signal. However, to be able to find contribution from at 
least three lines simultaneously we need higher amount of fre- 
quency channels than Planck HFI will have. A higher amount 
of frequency channels would certainly be possible for the next 
generation experiments like CMBPol even if they are based on 
already existing technology (Church 2002). 

There is the possibility that balloon and ground-based ex- 
periments will be able to check the level of enrichment of 
the universe by heavy elements even before Planck, observ- 
ing at I = 220 and I = 810, for instance. At these high mul- 
tipoles, effect is directly proportional to the optical depths in 
lines and the primordial CMB anisotropy. Hence using the data 
from tableland Fig|3(bottom), and the simple analytic rela- 
tion 6Ci ^ -2 Tjf, QP"'™''y ^ one can immediately give the ef- 
fect in a first approximation. The high signal-to-noise ratio of 
the primordial C;-s around the first three Doppler peaks will 
correspondingly give rise to strong signals in scattering, and 
might become accessible through the tremendous sensitivity 
promised from the forthcoming balloon and ground-based tele- 
scopes like Boomerang, ACT, APEX and SPT. 

CBI, VSA and BIMA interferometrs are studying success- 
fully CNB angular fluctuations at low frequencies 30 GHz, 
there are plans to continue observations on 45, 70 GHz, 90, 
100 and 150 GHZ. CI lines with lambda 609 and 370 micron 
coming from high redshifts 10-15 might contribute to the ob- 
served signal at low frequencies. 

Beginning of the end of 60-ties theorists are discussing 
early star formation due to isothermal perturbations, making 
possible production of heavy elements at redshifts well above 
100. This makes interesting the lines with much shorter wave- 
lengths, like 12.8// Ne II line and many others, which might be 
contributing to the Planck HFI spectral bands even from red- 
shifts ~ 100. 

5. Effect of Foregrounds 

Before WMAP, the ampUtude and spatial and frequency scaUng 
of foregrounds had not been firmly established, and its model- 
ing was merely in a preliminary phase. After their first mission 
year, WMAP's team have come up with a foreground model 
which is claimed to reproduce with a few percent accuracy the 
observed foreground emission (Bennett et al. 2003b). This and 
other future studies of foregrounds may provide a characteriza- 
tion of these contaminants such that their effect on our method 
can be minimized. We next proceed to estimate the impact of 
these contaminants on our method by the use of cuiTent fore- 
grounds models. 

We have adopted the middle-of-the-road model of Tegmark 
et al. 2000 (hereafter TOO). This model studies separately the 
contribution coming from different components, giving similar 
amplitudes for those which are also modelled by the WMAP's 
team, (Bennett et al. 2003b). In this model, we will consider the 
contribution of five foreground sources, namely synchrotron ra- 
diation, free-free emission, dust emission, tSZ effect (Sunyaev 
& Zel'dovich 1972) associated to filaments and superc lusters 
of galaxies, and Rayleigh scattering. The /-dependence of the 
power spectra was approximated by a power law for all fore- 



ground components, (except for tSZ, for which the model was 
slightly more sophisticated). The frequency dependence ob- 
served the physical mechanism behind each contaminant: sim- 
ple power laws were adopted for free-free and synchrotron, 
whereas a modified black body spectrum was used for dust. 
For tSZ, the frequency dependence of temperature anisotropics 
is well known in the non-relativistic regime. All details about 
this modeling can be found at TOO. We are neglecting the con- 
tribution from the SZ effect generated in resolved clusters of 
galaxies which can be removed from the map. We are also 
assuming that all resolved point sources are excised from the 
map, and that the contribution from the remaining unresolved 
point sources (cr^j in eq.(lO) of TOO) can be lowered down to 
roughly the noise level. This may require the presence of an 
external point-source catalogue. 

Fig|8^ shows the expected precision level when foreground 
contaminations can be neglected, and Fig 03 shows the effect 
of foregrounds in our differential method to detect the pres- 
ence of resonant species. Together with the power spectrum 
of our Standard ACDM model, (thick solid line at the top), 
we show the contribution from residuals of all foregrounds 
components under consideration, obtained after subtracting the 
HFI 100 GHz channel power spectra from the HFI 143 GHz 
one. The thin solid line corresponds to synchrotron emission, 
free-free emission is given by the dotted line; the dashed line 
gives the contribution of dust through vibrational transitions. 
The tSZ associated to filaments is shown by the triple-dot- 
dashed line, and the contribution from the combined noise of 
both channels is shown by the soUd line of intermediate thick- 
ness which crosses the plot from the bottom left to (almost) 
the top right corners. At these frequencies, the contribution of 
rotating dust is negligible, and the most limiting foreground 
is dust emission by means of vibrational transitions. Finally, 
the last frequency dependent contribution that we have con- 
sidered corresponds to Rayleigh scattering. As shown in Yu, 
Spergel & Ostriker 2001, hereafter YSO, Rayleigh scattering 
of CMB photons with neutral Hydrogen atoms introduce fre- 
quency dependent temperature fluctuations. However, provided 
that TRavieigh °^ , this proccss is only effective at high fre- 
quencies (v > 300GHz), causing deviations of measured power 
spectrum from the model C;'s of the order of a few percent. 
YSO characterized the Rayleigh scattering by introducing in 
the CMBFAST code a frequency dependent effective optical 
depth and a frequency integrated drag force exerted on the 
baryons. This last modification, which was neglected when 
considering resonant scattering, couples the evolution of the 
different temperature multipoles (Atj) for different frequen- 
cies. However, this coupling is exclusively due to the dipole 
term i, which, at the light of their results, is very simi- 
lar for all frequencies. Taking the same dipole term for every 
frequency allows us evolving a separate system of differential 
equation for each frequency. The frequency dependent changes 
in the angular power spectrum, (5C//C/(v)), obtained under this 
approximation are almost identical to those in YSO. For each 
channel pair under consideration, we computed the residual in 
the power spectrum difference due to these frequency depen- 
dence temperature anisotropics. We checked that this effect was 
always subdominant. In FiglSj^, the impact of Rayleigh scat- 
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tering is shown by a middle thick dashed line at the bottom- 
right corner of the panel, for the particular comparison of the 
lOOGHz and 143GHz channels. 

The addition of all contaminants (foregrounds + noise) is 
given by squares. Diamonds show the change in the power 
spectrum associated to a resonant transition placed at z = 25 
for optical depth t = 3 x 10 After comparing with Fig|8t 
we can appreciate how the presence of foregrounds affects our 
method in two different ways: a) it increases the minimum t 
to which the method is sensitive by about two orders of mag- 
nitude, b) it changes the range of multipoles to look at, due 
mainly to the large typical angular size of dust clouds. In the 
particular case of HFl 100GHz-143GHz channels, one should 
focus on the range I € [100,200] once the dust component is 
included rather than in / € [3,10] for the foreground-free case. 
Recurring again to the linear dependence of r versus the abun- 
dance, eq.@, we see that this model of foregrounds increases 
by about two orders of magnitude the minimum abundance on 
which future CMB missions will be able to put constraints. 

6. Conclusions 

In this paper, we have computed the effect of heavy elements 
on the Cosmic Microwave Background. We have selected those 
species which might leave a greater signature in the CMB 
power spectrum, and analyzed the possibility of constraining 
their abundances by the use of forthcoming CMB data. 

Under the approximation of negligible drag force, we have 
computed what are the changes in the CMB angular power 
spectra introduced by resonant species placed at redshifts 
which are dependent on the observing and resonant frequen- 
cies. The overall effect can be decomposed into a damping or 
suppression of the original CMB temperature fluctuations, and 
a generation of new anisotropics. For the optically thin limit to 
which we have restricted ourselves, damping dominates over 
generation of new anisotropics in the intermediate and high 
multipole range (I > 20). 

These changes in the CMB are of very small amplitude 
(0.01 - 0.3yu/r), and could be distinguished from the CMB 
component by means of their frequency dependence. Indeed, 
a comparison with respect to a reference channel containing 
only non-frequency dependent temperature fluctuations could 
be used to quantify the amount of angular power introduced by 
flie resonant species. However, this would only be feasible if ei- 
ther the other frequency-dependent components are negligible 
or characterized with extreme accuracy. The possible presence 
of foregrounds, galactic of extragalactic, whose amplitude and 
spectral behavior is still under characterization, is a serious as- 
pect to be taken into account. Other technical challenges, such 
as the calibration of different channels and the enormous sen- 
sitivity required, should be accessible from the next generation 
of detector technology. 

We have obtained limits for abundances of heavy elements 
when complete foreground removal is possible, but also have 
shown values when all foregrounds are present in the sky map. 
Our analyses have particularly focused on PLANCK HFI chan- 
nels, whose very low noise levels give extremely strong limits 
on abundances. It is easiest to observe the proposed effect with 



HFI at angular scales of 6* ~ 20° ox I ~ 10 because of the very 
low noise of the first three HFI channels at these multipoles. 
Using the 143 & 217 GHz channels of HFI, with the 100 GHz 
channel as reference, limits as low as 10"^ - 10"^ solar abun- 
dance were obtained for atoms and ions of the most important 
elements like carbon, nitrogen and oxygen in the redshift range 
[5,30]. At higher multipoles (/ > 200), we have shown that 
future balloon and ground-based experiments like ACT, APEX 
and SPT can put similar constraint on abundances, where the 
predicted signal is stronger because of the higher amplitude of 
the primordial CMB signal, and effect will be directly proporti- 
nal to the optical depths in lines. The presence of foregrounds 
makes all these limits about a factor of 10^ worse, but that must 
be considered as the most pessimistic scenario. 

In this paper we have considered scattering of CMB pho- 
tons by low density gas everywhere in the universe, where op- 
tical depths in the relevant atomic or molecular lines are very 
low and resonant scattering effect dominates over line emis- 
sion. The effect connected with scattering from atomic and 
ionic fine-structure transitions depends on the amount of of 
atoms/ions in the ground state, which is determined by the ex- 
citation temperature. In Appendix B we showed that even up to 
very high over-density of the plasma, the change in excitation 
temperature from the CMB temperature is negligible, and our 
simple approximation is valid. The clustering of heavy atoms 
and molecules in clouds should change the predictions on the 
modifications of the CMB power spectrum at smaller angular 
scales, and will be discussed in a subsequent work. 
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Appendix A: 

METHOD OF COMPUTATION: 

Our approach to model the coupling between these heavy 
species and CMB photons will be similar to that outlined in 
Zaldarriaga & Loeb (2002, ZL02). In their paper they dis- 
cussed modifications of CMB spectrum by presence of pri- 
mordial Lithium atoms. Here we extend their analysis to other 



atoms and ions. We also discuss the necessary changes that are 
enforced while extending this method to all resonant species. 

As a first step, we describe the modifications introduced in 
the CMBFast code in order to compute the effect of resonant 
transitions. We consider a given resonant transition / of a given 
species X, with a resonant frequency vx, ■ For a fixed observing 
frequency Vohs, the redshift at which that species interacts with 
the CMB is 1 -H z = vxjvohs, and its opacity can be written, in 
general, as fx, - Tx^Pir]), with a normalized profile function 

dri'Virj') - 1. We shall model this profile with a gaussian: 



exp 



Icrl 



(A.l) 



where Tx, is the optical depth for the specific transition, rjx; 
is the conformal time corresponding to the redshift 1 + zxj - 
Vx, I Vobs where scattering takes place, and ctx, is the width of the 
gaussian. For a fixed transition, this width should be given by 
the thermal broadening of the line. For the sake of simplicity, 
we shall take crxj ~ O.OIt/x, . 

Once the line optical depth has been characterized, (as 
shown in Section 2), the new opacity term with a gaussian pro- 
file is added to the standard Thompson opacity. 



T - aneCTT + 



Tx,' 



4- 



27TO-1 



(A.2) 



and this in turn is used to compute the visibility function 
T(7;), defined as 



T(7,) = f(77)e 



(A.3) 



This scattering introduces a frequency dependent term in 
the evolution equation for the photon distribution function, 
which results in a frequency dependence drag-force. However, 
based on the same arguments of ZL02 , we can safely ignore 
this drag-force as long as the characteristic time of drag 
exerted by these species is far larger than the Hubble time, 
which is indeed the case due to the low optical depths under 
consideration. These are the modifications we introduced in 
the CMBFAST code in order to compute the CMB power 
spectrum under the presence of scattering by atoms and ions. 

ANALYTICAL FORM OF SC, 's: 

Following eq.0 for the A;-mode of the photon temperature fluc- 
tuation, 

rno 
Jo 

X [r(T])(ATQ-ifiVb) + ^ + ip-ikfii^l (A.4) 

we proceed now to characterize the change in the radiation 
power spectrum. In ea.( IA.4> . the terms in the angle yu = k ■ n 
can be eliminated after integrating by parts and neglecting the 
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contribution of boundary terms, (Seljak & Zaldarriaga 1996). the profile peaks at rj - rjXj and that and 7;+ are such that 

Virf) ^ for T] < rj- and rj > Hereafter, tjx, will be referred 
to as transition epoch or line epoch. If we add this new term to 
(A. 5) the opacity, and assume that t^, <k 1, then we can expand the 
term SAjj in a power series of ■ In this case, we obtain: 



This gives: 
rm 

At = dT]e''"'^'^-'^''^S(k,T]), 
Jo 

where the source term S (k, rf) is defined as 

S(k, if) = (0 + «a) + T (a™ + «A + y ) + T (y ) ■ (A.6) 

From this source term, the angular power spectrum can be 
expressed as, (e.g. Seljak & Zaldarriaga, 1996): 



^At-,, = tx, ■ ^ - j dri ji [k{r]Q - ?])] ^(77) 



Ci = (4: 



= (47r)2 J 



dk k^P^(k) 



rm 2 

I dT]S (k, rj)ji [k{T]u - T])] 
Jo 



dk k^P^ik) \Atj\ 



(A.7) 



(0 + i/r) + re ^At-o + i/r + y j + 

r d7]ji[k(Tjo-Tj)]'P(T])e-' 
Jo 

X (4„, 2t)- 



where ji is the spherical Bessel function of order I and P^(k) 
is the initial scalar perturbation power spectrum. If at a given 
frequency the CMB interacts through the resonant transition / 
of a species X, the measured power spectrum will differ from 
the reference one by an amount: 

6Ci = Cf' — Ci 



(47rf J dk k^P^ik) [ dri S^'-iKTi)], [^(770 - 77)] 
- dr] S(k,7])ji[k{T]o-7])] J 



+ Tv 



^ J^'* d?J j,[k{T]o-T])]je(Tj) 

(0 + 1^) + te"^ ^Ato + (ff + + ^ {'^ + ■'"^) 
r'l+ 

+ di] ji[k{7]o-7])]e-^ X 
Jo 

-(Aro + <A+y)wn^) 



= (Any 



- At/ 



71)2 J k^p^(^k)[2 At,! + 5At,i]6Atj (A.8) = tx, ■ Di + t% ■ Di + 0{tI). 



X |(^2(^)+y[(;7)(-^ 



2f 



+ 0{tI) 



(4: 



with 5Ar,/ = X"" t/77 77) - 5 (k, Tf)) ji {kirio - rf}], and 

where the term S^'{k,rf) refers to the sources including the 
species X. Note that the cross term 2 Ajj ■ 6Aj i only arises af- 
ter computing the difference of the power spectra, i.e., it is not 
present if one computes the power spectrum of the difference of 
two maps obtained at different frequencies. This cross term is 
precisely the responsible of having dCi Unear in Tx,, and hence 
also linear in the abundance of the species X, . This term also 
accounts for the correlation existing between the temperature 
fluctuations generated during recombination and those gener- 
ated during the scattering with the species X. This non-zero 
correlation is essentially due to those k modes corresponding to 
lengths bigger than the distance separating the two events, i.e., 
recombination and scattering with X, (Hernandez-Monteagudo 
& Sunyaev, in preparation). Let us now model the differential 
opacity due to this transition as tx, = Tx^Virf), where !P(77) is a 
profile function of area unity, drj' Virj') - 1 . We write the 
optical depth as TtoArf) - T{ri)+TxXv), where the last term refers 
to the optical depth due to the transition / of the X species. It 
can be expressed as Tx^iri) - tx, -^(77) — tx, J^" dri' P(ri'), with 
^(77) the area function of the profile 'P(77). We will assume that 



(A.9) 



It is worth to note that Di contains both the suppression of 
intrinsic CMB anisotropics, and the generation of new fluctua- 
tions at line epoch. This expansion allows us writing 6C1 as a 
power series of tx, as well, 

SC, = Tx, ■ (47r)2 J dk k^P^k) [Di 2Atj] 

(47r)2 j dk k^P^ik) [D2 2Atj + Dj] 



+ Ty 



(A. 10) 



The result of this expansion is displayed in Fig lA.ll the 
change in the power spectrum induced by resonant transitions 
is computed for different redshifts, 25 (top panel) and 500 
(bottom panel), and different amplitudes of the optical depth 
(solid lines correspond to t = 1.5 x 10"^ and dot-dashed lines 
to T = 0.15). Diamonds give the linear approximation in r, 
whereas triangles show the quadratic one. Both match the ex- 
act 6Cfs fairly well for the low t cases. Therefore, by means of 
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eq.Q, we can establish a linear relation between 6Ci and the 
abundance of the species. 

In Fig |A.2l diamonds show the absolute value of C\ versus 
I. As shown above, this term is the sum of two integrals. The 
first one is the suppression induced by the e '"'^^ term, (thick 
solid line in the figure), whereas the second contain the newly 
generated anisotropics, (thick dashed line). The latter term has 
a monopole (Aq + tfr, thin dashed line) and a velocity (thin dot- 
ted line) contribution. We are plotting absolute values of all 
terms. As pointed out by ZL02, the monopole term decreases 
with cosmic time and hence the velocity term becomes the most 
important source of generation of new anisotropics. However, 
this generation takes place at the transition epoch, and hence is 
shifted towards the low multipole range. It is easy to show that 
the maximum multipole below which generation becomes sig- 
nificant is roughly given by Zx, ~ 27r (rjo - rjx^ /rjx^- We show 
that the suppression term is dominant for high multipoles, and 
only at very low multipoles suppression and generation tend to 
cancel each other. Only if higher orders in the power expansion 
are relevant, (i.e., if tx, ~ 1) the newly generated anisotropics 
become important. 

This cancellation of suppression and generation terms at 
low multipoles can be better understood when coming back 
onto ea lA.4l If in this equation we substitute t by r + tx, , we 
find that the change in the temperature modes due to tx, can be 
written as: 

rno 

6 At = (e-^"' - 1) drj e'*f (''-'»)te-" (Aro - /yuv/,) 

Jo 

r'lo 

+ e-'': di] e*^("-'">>fx,e"^ (Aj-o - ilivb) , (A.ll) 
Jo 

where we have applied the definition of the visibility func- 
tion, T(77) - f(j]) e"'"*'''. However, recalling that the visibility 
function gives the probability of a photon being emitted at a 
given 7], and taking fx, = Tx^Soil-riXi), the last equation merely 
implies that 

6 At ^ Tx, [A7-o(77A',) - ATo(7]rec) ' (^'fe('7Z,) - Vb(T]rec))] , (A. 12) 

where we have taken into account that tx, is much smaller 
than unity, and approximated the exponentials to unity, as 
we shall focus on the very low k range. That is, the change 
in the temperature mode AT{k,fj.,ri) reflects the difference of 
the monopole (Aj-q) and velocity terms when evaluated at 
recombination and at the transition epochs. We can estimate 
how this difference projects onto multipole space if we restrict 
ourselves to the low multipole (large angle) range. For the 
scattering redshifts under consideration (zx, > a few), we can 
safely neglect the term due to the Integrated Sachs-Wolfe 
effect. Then from ea.( IA.4> . one finds that the time depen- 
dence of the monopole term AToirj) can be approximated as 
A7-o(?7) - ATo(T]rec) jo (k [j] - J]rec]), whcrc jo(x) is the spherical 
Bessel function of order zero. From this dependence, it is easy 
to see that at low multipoles, i.e., at small enough k's, A7-o('7m) 
and AToirjriy.) will be roughly equal, and thus their difference 
in the equation above will tend to vanish. This can also be seen 
in Fig lA.2l where the contribution of the monopole term to C\ 
(thin dashed line) has roughly the same amplitude at redshifts 



500 and 25 for Z < 10. This behaviour is the responsible for 
the cancellation of the 5C/'s at low /'s. In this situation, the 
difference of the velocity terms will be of particular relevance. 
The evolution of velocities can be approximated, after inte- 
grating in ^l, as Vh(ri) - Vb(j]rec)D{7])/D{j]rec)j\ (k [t] - Tlrec]), 

with D(ri) the linear growth factor . The growth of velocities 
will assure an increase of the Doppler-induced anisotropy 
power. Finally, due to the fact that ji(x) is maximum at 
X ~ I, then we have that the Doppler term will predominantly 
contribute for '' k ~ In/rix,, which corresponds to multipoles 
Zx, ~ 27r (?7() - ?7i ) /t/x,- Again, this can be checked by looking 
at the velocity term (dotted line) in Fig lA.2l the power is 
projected at lower multipoles at later epochs, scales at which 
the amplitude grows with conformal time as D. 

Therefore, the factors determining the cross-over of the 
5C/'s from (Doppler-induced) positive values at small / to 
(absorption-induced) negative values at large I are two: a) the 
constant amplitude of the monopole at large scales, and b) the 
growth of peculiar velocities with cosmic time . The angular 
scale at which such crossing takes place is roughly determined 
by the time at scattering epoch, t/x,, /x, ~ 27r(77o - ??Ai) /'7a', • 



Appendix B: 

LINE EMISSIVITY: 

The volume emissivity for line emission due to collisional ex- 
citation can be written as 

,. hy 

X"" = — <«,■ «c) r/« (^(v'), (B.l) 

expressed in erg cm^^s^^Hz^^str^^ in the rest frame of the 
emitter. Here n,- is the number density of the atom or ion under 
consideration, and ric is the number density of the colliding par- 
ticle, and (n, tic) means averaging over the volume specified by 
multipole /, corresponding to the angular size of observation, 
and the thickness of the slice along the line-of-sight defined by 
the frequency resolution of the experiment, hv is the energy of 
the emitted photon, y/„ is the collision rate-coefficient between 
between the two levels, and (^(v') is the fine-profile function, 
which we can approximate by a 5-function because experimen- 
tal bandwidth is much larger than velocity or thermal broaden- 
ing, 

^(v') = S(v' - v) = — 6(z' - z) , (B.2) 

V 

where z corresponds to the redshift under study. Fan. ( IB. it 
is correct for the low density gas (« <K rtcri,) that we are con- 
cerned with in this paper. This emissivity produces a line flux, 
written as following in an expanding universe (neglecting ab- 
sorption because Ty <s; 1) 

- j JYTz)^ dl, (B.3) 

We are assuming that zyx, >> tj^^c throughout the paper. However, 
strictly speaking, one should consider the difference of conformal 
times of recombination and resonant scattering, i.e., k ~ 2n/(T]xj -J]rec)- 
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where the integration is along the line-of-sight. This creates 
a distortion in the thermal spectrum of the background CMB 
photons 



AT e' -I Fy hv 
't~ ~ «~ R~ ' where x = — — 



(B.4) 



with Ty as the mean temperature of CMB photons, and By 
being the Planck function. 



EFFECT OF OVER-DENSITY: 

In this subchapter, we consider the impact of collisions in our 
study. Collisions with electrons (or neutral atoms in the case of 
neutral gas) cause emission from the same atomic and ionic 
fine-structure lines and hence change the observed <5C/-s by 
producing an additional and independent signal on the top of 
the scattering signal. 

In addition, collisions are decreasing the population of the 
ground level in the transitions under consideration, dimishing 
the optical depth due to scattering in the line. In what follows, 
we compute the deviation of the equilibrium excitation temper- 
ature of a two-level system from the background cosmic mi- 
crowave radiation temperature as a function of the over-density 
of the enriched region, we estimate the required over-density 
for causing a given amount (30% or 50%) decrease in the opti- 
cal depth. 

We consider a two-level system (like CII or NIII ion) for 
simplicity, whose population ratio between the two levels is 
determined by the excitation temperature, Tex, defined as 



(n„/n/) = {gulgi)&y.p{-hvlkTEx\ 



(B.5) 



where n„ and «/ are the relative fraction of atoms/ions in 
upper or lower levels, respectively, and g„ and gi are the statis- 
tical weights of each level. The equilibrium population ratio is 
obtained by solving the statistical balance equation: 

«„ (A„, + B,iJy + n,y„i) = n, I^B.iJy + i^n,j,,e-'"'"^A (B.6) 

gi I 

Here A^i and B„; are the Einstein coefficients for sponta- 
neous and stimulated emission, and Jy is the background radi- 
ation field of CMB photons, Jy = (2hv^ / c^)[exp{hv / kTcMB) - 
1] TcMB is the radiation temperature at redshift z, defined 
by TcMB = To(l + z), with Tq = 2.726 K. is the elec- 
tron density inside the object at redshift z , connected to its 
over-density 6 by the following: rieiz) = he{z) (1 +6), where 
heiz) - 2.18 X 10^^(1 + z)^ cmT^ . jui is the collisional de- 
excitation rate (in cm^ s"'), and Tk is the kinetic temperature 
of the electrons, which we assume to be roughly constant at 
10"* K. Then from ean.( IB.5t we can express the over-density as 
a function of excitation temperature: 



A,i [(1 + r)e-'"'/*^^'' -r] 



neiz) Jul [e 



-hv/kTK 



-hvjkTE. 



(B.7) 



where T = [exp{liv / kTcMsiz)) - 1] and we have used 
Bui - (c^ /2hv^)Aui. Eqn. JB.7t is exact, which can be further 
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Fig. B.l. Over-density and actual electron number density re- 
quired to cause significant deviation of the excitation temper- 
ature from the background cosmic microwave temperature, as 
function of redshift. Shown here are the cases for CII and NIII 
fine-structure singlets. The thick lines are for CII ion,and the 
thin ones are for NIII ions. The pair of solid and dashed lines 
correspond to 30% and 50% decrease in the optical depth for 
the same column density, respectively 

simplified in the limit of low redshifts, when the line is in the 
Wien part of the CMB spectrum {Tcmb{z) < Tex ^ hv) 



1 = 



«.(z)r<4i-;^-^"™'] 



(B.8) 



At rig < A,ii/y,ii collisions do not influence the population 
of the ground level because radiative transitions are faster than 
excitations due to collisions. If density is higher, population 
of the ground level decreases and optical depth due to reso- 
nance scattering drops correspondingly, diminishing the effect 
we have discussed in this paper 

Obviously, a very large over-density is needed in order to 
produce any significant amount of deviation from the CMB 
equilibrium, as can be seen from Fig lB.ll for the case of CII 
158/z and NIII 57// transitions. These two transitions are picked 
as they both arise from two-level fine-structure splitting, but 
results are similar for all other relevant atoms and ions. In the 
bottom panel of Fig lB.ll we see that the electron density re- 
quired to cause a given amount of depopulation of the ground 
level is constant for low redshifts, and this density is less than 
the critical density for that transition (e.g., with CII 158jU line, 
we have y,,/ = 4.6 x 10"^ cm-' s"' for collision with electrons at 
lO'^K, and the ratio Auiljui - 50 cm This simple estimate 
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gives 5-10 times higher critical density than the more exact 
approach shown in Fig lB.U . Hence, we conclude that the ef- 
fectiveness of the resonant scattering process decreases in the 
most over-dense regions of the universe. Atoms and ions of 
heavy elements everywhere else including less over-dense ha- 
los and clouds participate in the process and contribute into the 
angular fluctuations of the CMB we discussed in this paper. 
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Fig. 3. Top: A schematic diagram for the abundance history of the universe, (a) The upper panel shows the global pattern of 
metal abundances, showing the two major epochs of enrichment of the IGM: first, during the peak of activity of massive Pop 
III stars around redshift 15 - 25 (shaded area), and second, during the peak of galaxy formation around redshift 3-5 when 
global rate of star formation reaches maximum, (b) Middle panel shows the relative fraction of three major atomic and ionic 
species: OI (or HI), CII and OIII, normaUzed so that the total abundance of all ions of a given element is close to unity at any 
redshift. The OI abundance closely follows the neutral hydrogen fraction of the universe because their almost similar ionization 
potential (see discussion in text). The redshift scale is chosen as log(z + 5) to emphasize the redshift region 10 - 30 of interest 
to this paper At very low redshifts (z < 0.7) the IGM gets heated to very high temperatures (T ~ 10^ - 10^ K) causing even 
higher ionized species to exist, e.g. OIII OVI. (c) Lower panel shows another variation of relative ion fraction, where Pop III 
stars ionize all the oxygen around redshift 15, so that OIII have higher abundance and correspondingly OI and CII have lower 
abundance. Bottom: Frequency dependence of optical depth in atomic and ionic fine-structure lines, in accordance with the 
abundance history sketched above. Shown here are the contributions from the four most important lines: neutral oxygen 63jU, 
doubly ionized oxygen 88;U, singly ionized nitrogen 205/i and singly ionized carbon 158jU, and their total contribution for each 
histories. History of Nil ion is taken as similar to that of CII ion. 
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Fig. 4. Frequency dependence of temperature anisotropy generated by scattering from fine- structure lines. We take two ionization 
histories and the high abundance case from figure and show the contributions of different lines in different frequency range, for 
a fixed multipole I - 10. Four different redshift ranges are marked for each species to emphasize the epochs where the dominant 
contribution from each lines are coming. The sensitivity limits of Planck HFl channels are marked by the crosses. The HFl limits 
have been improved by a factor of VAZ by averaging over in the multipole range I - 7-16 . The sensitivity limits are for 3 standard 

deviation detection, and the y-errorbars correspond to Icr error in 3cr detection (we recall that cr — 3 -yJcQiprobe) '^Ci(refy where 
the reference channel is fixed at HFl lOOGHz). The x-errorbars corresponds to the wide bandwidth (~25%) of Planck HFl 
channels. 
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Fig. 6. Temperature anisotropy generated by atomic and ionic fine structure lines at different frequencies in the case of three 
angular scales of observation. The figures on left are for Case I of Fig. 3, which is for reionization by stars unable to ionize OIL 
Figures on right are for Case II, which is for reionization by hot stars and quasars and hence have high OIII fraction. Each figure 
gives the total contribution for both history A and history B. However, the separate contributions from difl'erent Unes are given 
for history A only. We present results for three angular scales: Z = 10 (top), I - 220 (first Doppler peak, in middle), and Z = 810 
(third Doppler peak, at bottom). The scattering-generated 6Ci-s are proportional to the intrinsic C/-s for higher multipoles, but 
magnitude is sUghtly less at the position of third peak because of the reduced power in primary spectrum from / = 220 to / = 810. 





545 353 217 143 545 353 217 143 

v„hs (GHz) y„4,, (GHz) 



Fig. 7. Late ionization and enrichment history of the universe. Top, left: Sketch for the abundance history with late reionization, 
where first metals are produced around z ~ 10 - 12. Top, right: Temperature anisotropy for this abundance history generated at 
Z = 10. Also shown are the Planck HFI sensitivity, where the noise levels have again been improved by VaI in the multipole 
range 10 - 20. Bottom left and right: Temperature anisotropy for / = 220 and / = 810, respectively. 
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Fig. 8. {a) Example of optical depth due to a resonant transition which can be detected by the HFI 100 GHz and HFI 143 GHz 
PLANCK channels in the absence of foregrounds. The upper thick solid line gives the reference model CMB power spectrum. 
Diamonds show the expected difference in the power spectrum from both channels due to the resonant scattering placed at 
z = 25 with T = 1 X 10"^. {b) Presence of foreground contamination after subtracting HFI lOOGHz channel from HFI 143GHz 
channel. The reference ACDM CMB power spectrum is shown in thick solid Une, whereas the instrumental noise after the 
map subtraction is shown in solid intermediate-thickness line. All thin lines refer to the foreground model as quoted as in TOO: 
vibrating dust emission is shown by a dashed line, whereas free-free and synchrotron are given by a dotted and a sohd thin 
Une respectively. tSZ effect associated to filaments and superclusters is given by the triple-dot-dashed line. Rotating dust gives 
negligible contribution at these frequencies. Rayleigh scattering introduces some frequency-dependent variations in the C;'s, but 
well below the noise level (thick dashed line at the bottom right corner). Diamonds show the amplitude of the change in the CMB 
power spectra induced by resonant species placed at z = 25 and r = 3 X 10"^. 
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Fig. A.l. This figure shows the vaHdity of the linear (diamonds) and quadratic (triangles) approximations on t when trying to 
describe the changes induced in the power spectrum by a resonant transition (5C;'s). We see that both suffice to accurately match 
the theoretical 5C;'s for redshifts 25 (top panel) and 500 (bottom panel) and for t - 1 .5 x 10"^, (solid lines). However, when t is 
closer to one (dot-dashed lines), the quadratic approximation performs remarkably better than the linear one. 
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Fig. A.2. Study of the different contributions to Ci in eq. (ll2> . Diamonds give Ci versus multipole. We are plotting only absolute 
values. In units of r^,, we plot the damping term in solid line. The term responsible for the generation of new anisotropics is 
plotted in thick dashed line. Its two components, monopole ((Aq + ij/)) and velocity, are also displayed in thin dashed and thin 
dotted lines respectively. As conformal time goes by, the relative weight of the velocity term with respect the monopole term 
increases, and the sum of both at low multipoles approaches the amplitude of the original CMB anisotropics, (solid line). Both 
generation and damping tend to cancel each other at a multipole I - /x,, dependent on the epoch of interaction rjx^. 



